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Abstract 

Strong-field ionization generally produces ions in a superposition of ionic eigenstates. This su- 
perposition is generally not fully coherent und must be described in terms of a density matrix. A 
recent experiment [E. Goulielmakis et al, Nature 466, 739 (2010)] employed attosecond transient 
absorption spectroscopy to determine the density matrix of strong-field-generated Kr"*" ions. The 
experimentally observed degree of coherence of the strong-field-generated Kr+ ions is well repro- 
duced by a recently developed multichannel strong-field-ionization theory. But there is significant 
disagreement between experiment and theory with respect to the degree of alignment of the Kr"*" 
ions. In the present paper, the theory underlying attosecond transient absorption spectroscopy of 
strong-field-generated ions is developed. The theory is formulated in such a way that the nonper- 
turbative nature of the strong-field-ionization process is systematically taken into account. The 
impact of attosecond pulse propagation effects on the interpretation of experimental data is inves- 
tigated both analytically and numerically. It is shown that attosecond pulse propagation effects 
cannot explain why the experimentally determined degree of alignment of strong-field-generated 
Kr"^ ions is much smaller than predicted by existing theory. 

PACS numbers: 32.80.Rni, 82.53.Kp, 31.15.A-, 42.65.Re 
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I. INTRODUCTION 



The interaction of matter with a hght pulse sets electrons in motion. The associated 
dynamics frequently occur on time scales comparable to or smaller than a femtosecond. The 
time-domain observation of such fast phenomena was impossible until a few years ago, when 
the first pump-probe measurements with attosecond extreme-ultraviolet (EUV) pulses were 
demonsuated fl. Siuce the. several types of attosecoud spectroscopy have been estabhshed 
or proposed theoretically Most of these approaches employ an intense near-infrared 
(NIR) laser pulse that either induces or probes dynamics of interest. Combined with an EUV 
pulse, an attosecond pump-probe measurement is performed by varying the delay between 
the pulses and observing the results of their interaction with a system under study. One 
of the most prominent examples is attosecond streaking js-S], where photoelectron spectra 
are measured for a set of delays. Attosecond streaking has proved to be a very powerful 
technique for measuring electron dynamics triggered by an attosecond pulse |6|, but it is 
less suited to probing strong-field dynamics with an EUV pulse. Historically, the first time- 
resolved measurements of electron dynamics induced by a strong light field were performed 
with the aid of attosecond tunneling spectroscopy [3], where ions, rather than electrons, 
were measured for different delays between the EUV and NIR pulses. Very recently, another 
promising technique — attosecond transient absorption spectroscopy — has been established: 
by measuring the transmission of an EUV pulse through a gas ionized by an NIR pulse, the 
motion of ionic valence-shell electrons has been observed with attosecond resolution 8|. 

The earliest implementation of transient absorption spectroscopy for the study of ultra- 
fast dynamics involved detecting the change in absorbance of a narrowband probe pulse as 
a function of pump-probe time delay j9-12|. In experiments that elucidate ultrafast molec- 
ular dynamics, tuning the central wavelength of the probe pulse tracks the motion of the 
nuclear wave packet into and out of different regions of the excited-state potential energy 
surface. Advances in laser pump-synchrotron probe techniques have allowed time-resolved 



x-ray absorption spectra to be collected in this fashion 13 



Mathies and Shank introduced a variant of transient absorption spectroscopy in which 
the photoexcitation of a sample by an optical pump pulse is followed by probing with a 

o|4 



spectrally broadband, few-cycle visible pulse 



20N23|. The probe pulse that is transmitted 



through the sample is spectrally dispersed by a spectrometer. In this method, the time 
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resolution is mainly determined by the duration of the pump and probe pulses, and it is in- 
dependent of the spectral resolution. Hence, compared to earlier approaches that employed 
narrowband probe pulses, spectrally resolved transient absorption spectroscopy with broad- 
band probe pulses offers the obvious advantages of high time resolution and high spectral 
resolution. In the visible, state-of-the-art experiments on samples of biological relevance have 
been done with sub-5-fs pump and probe pulses [2J| . This technique has also been applied 



to femtosecond studies of novel materials 



25 



26|. For an up-to-date summary of the devel- 



opment of the optical transient absorption technique, see Ref. [27[. Recently, femtosecond 
transient absorption spectroscopy was extended to shorter probe wavelengths by employing 
laser-produced bremsstrahlung x-rays 
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as femtosecond synchrotron x-rays 



and EUV high-order harmonics |29|-|32| . as well 



33|. 



Transient absorption spectroscopy provided direct evidence for hole alignment in strong- 
field-generated atomic ions 15|. The ion alignment dynamics driven by electron- ion col- 
lisions in a strong- field-generated plasma were studied in Ref. 3J]- However, it was not 
clear at the time whether strong-field ionization leads to the formation of a coherent super- 
position of ionic eigenstates (in cases where not only the ionic ground state is populated). 
In other words, it was not clear whether strong-field-generated ions undergo any coherent 
intraatomic dynamics. Theory indicated that it is generally not possible to describe strong- 
field-generated ions in terms of a perfectly coherent wave packet; a description in terms of a 
density matrix is required 35|, |36[ . Using transient absorption spectroscopy in combination 



with theory, the dia gon al elements of the density matrix of strong-field-generated atomic ions 



were measured 



16 



29[. The diagonal elements of the density matrix (represented in the 



ionic eigenbasis) are the populations of the various ionic eigenstates, generated by a strong 
NIR field. The measurement of the entire ion density matrix — including the coherences, 
i.e., the off-diagonal elements of the ion density matrix — was made possible by attosecond 
transient absorption spectroscopy Q]. 

In the present paper, we develop a theoretical description of attosecond transient absorp- 
tion spectroscopy of strong-field-generated ions. Our approach allows us to treat the pump 
step, i.e., the interaction with a strong NIR field, in a completely nonperturbative fashion. 
We assume throughout that the attosecond EUV probe pulse has no temporal overlap with 
the strong NIR pump pulse. In Sec. [Tll we derive the basic expressions underlying the 
theory. The analysis of the experiment of Ref. 8| becomes particularly transparent if it is 
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assumed that the attosecond probe pulse remains sufficiently short as it propagates through 
the NIR-modified target medium. This short-pulse approximation is discussed in Sec. Illli 



ow us to assess the validity of 
The distortion of the NIR 



Numerical EUV pulse propagation calculations in Sec. [IV] al 
the short-pulse approximation under the conditions of Ref. 
pump pulse by the target medium is not considered. Section |V] concludes. Atomic units are 
used throughout, unless otherwise noted. 



II. GENERAL CONSIDERATIONS 



A. Atomic response 

We consider a semiclassical description of the interaction of the pump and probe pulses 
with the atoms in the target gas. Both pulses are approximated by transverse, infinite plane 
waves propagating along the x axis, and are assumed to be linearly polarized along the z 
axis. As mentioned in the introduction, the NIR pulse is assumed to remain unmodified as 
it propagates through the gas. Thus, neglecting diffraction, the NIR electric field may be 
written as i^NiRl^^L — x/c), where is the time measured in the laboratory frame and c is 
the vacuum speed of light. For an atom at position x, it is convenient to introduce the local 
time t = ti — x/c. Hence, in the electric dipole approximation, the Hamiltonian for an atom 
at position x reads 

H{t) = Ho-Eo- £mR{t)Z - Sevy{x, t + x/c)Z. (1) 

Here, Hq is the unperturbed atomic Hamiltonian, Eq is the atomic ground-state energy, Z 
is the z component of the electric dipole operator, and £^EUv(a;, ^l) is the EUV electric field. 
In order to calculate the EUV-induced polarization response of the atoms, we need to solve 
the time-dependent Schrodinger equation 

z^^\^,t) = Hm,t). (2) 

Let us assume we have solved the x- independent, NIR-only problem 

i^UmRit, -oo) = ji^o - Eq- £mR{t)Z^ UmRit, -oo). (3) 

A suitable initial condition for the time evolution operator is 

f^NiR(^, —oo) — !■ exp |— ?(_^o — -E'o)^! as t — )■ — oo. (4) 
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Therefore, in the absence of the probe pulse, a solution to Eq. ([2]) is 

\^mR,t) = UmRit, -oo)|^o)- (5) 

Here, |\l/o) is the initial state of the atom, assumed to be the ground state. 

In order to take into consideration the effect of the probe pulse, we make the ansatz 

|*,t) = |*NIR,t) + |*,t)^'^ + ..., (6) 

where t)^^^ is a correction that is of first order with respect to £E\jv{x,t + x/c). It follows 
that 

i-\^,t)^'^ = [Ho - Eo - £mR{t)Z^ \^,t)^''^ - SEvv{x,t + x/c)Z\^^i^,t). (7) 
This equation can be integrated analytically. The result reads 

\^,t)'-^^=t[ dt'UmR{t,t')ZSEvy{x,t' + x/c)\^mR,t'), (8) 

J — oo 

where UmRit,t') = UmR{t, -oo)f/J,jj^(t', -oo). 

Let us now calculate the polarization along the z axis: 

P{x,t + x/c) = nAT{^,t\Z\<i/,t) (9) 
= Pncix, t + x/c) + P(^) {x, t + x/c) + .... 

In this expression, uat is the atomic number density, 

Phg{x, t + x/c) = nAT(^NIR, tl^l^NIR, t) (lO) 

describes harmonic generation driven by the NIR pulse (no EUV pulse present), and 

P^'\x,t + x/c) = nAT{^mR,t\Z\^,t)^'^ + c.c. (11) 

is the polarization correction to first order with respect to SEvy{x,t + x/c). 
Using Eq. ([8]), the first-order polarization correction may be written as 



P^^\x,t + x/c) = iriAT I dt' \SEvy{x,t' + x/c) 

X (*NIR, t\ZUmRit, t')^|*NIR, t') } + C.C. (12) 



This result, which is vahd for arbitrary pump intensities, generahzes Eq. (17b) in Ref. 22 1 



and shows that Eq. (21) in that paper is not quite correct 45|]. Note that the integrand in 
Eq. (IT^ depends exphcitly on t via (\l/NiR,t| and UNm(t,t'). Therefore, even if the EUV 
pulse may, effectively, be approximated by a delta function centered at, say, t' = t-EVY, the 
polarization induced by the EUV pulse contains information not only on atomic properties 
at the instant of the EUV pulse, but, in principle, also on atomic properties after the EUV 
pulse (for t > tEuv)- As we will see, this causes no difficulty for the situation considered in 
this paper. 

We assume that the EUV probe pulse comes after the NIR pump pulse, so that the NIR 
pulse does not affect the electronic states reached via EUV photoabsorption. Hence, we 
have, for the time evolution operator in Eq. f lT^ . 

f/NiR(t,t') = e-'(^«-^«)(*-*'\ (13) 

Let |/) denote an eigenstate of the unperturbed atomic Hamiltonian Hq with eigenenergy 
Ej. We may then expand the NIR-only state vector, after the NIR pulse, as follows: 

The expansion coefficients aj are time-independent. Thus, Eq. (IT^ goes over into 

1,1' F 

X r dt'£Euv(x,t' + a;/c)e-^(^^"^^)(*-*')+*(^^-^^')*'+c.c., (15) 

J — oo 



where \F) is an eigenstate of Hq with eigenenergy Ep, which is, in general, assumed to be 
complex. More precisely, the imaginary part of Ep is either zero or negative. 

B. EUV pulse propagation 

We assume that the high harmonics generated by the NIR pulse [Eq. (fTOj) ] do not overlap 
with the spectral range of the EUV pulse. Hence, the propagation of the EUV electric field 
through the medium can be described by the following scalar wave equation: 

/ 92 1 92 \ ^ , , 47r 92 



With the Fourier representations 

SEVY{x,tL) = — |W(x,a;)e^^'^(*^^^-/'=) + c.c.|, (17) 

p('\x,tL) = {P(^)(^,^)e-^"^*^-^/^) + c.c.} , (18) 

it follows from Eq. f|T6|) that 

In the situation considered, the spatial derivative of the electric field amplitude at a given 
u changes slowly over a wavelength 2ttc/u. This allows us to neglect in Eq. (fT9|) the second 
derivative with respect to x. Thus, the differential equation governing the spatial evolution 
of SEvy{x, ^) reads 

—SEvy{x,uj) = 27rz^P«(x,a;). (20) 
Equation (120|) can be integrated analytically if P^^\x,uj) is proportional to £evy{x,oj)- 

P^'\x,u;)=x^'\x,u)SEvy{x,u). (21) 

Here, x^^Hx^co) is the linear susceptibility. If Eq. (12T]) is valid, it follows from Eq. (120|) that 

£euy{x,uj) = £EUY{xo,uj)exp i^27ri'^ j dx' x'^'^\x' , uj)^ . (22) 

In the attosecond transient absorption experiment described in Ref. js], the EUV radiation 
transmitted through the sample (length L) was spectrally dispersed and analyzed. The 
detected signal is then proportional to 

\SeML,u:)\' = |£~Euv(0,a;)|V^-^/o^'^^i-[>^'^'(^'-)l, (23) 

which, for a homogeneous target medium, goes over into 

\£evy{L,co)\' = |4uv(0,a;)pe-^-^^^-[>^'^^(-)l. (24) 



Equation (!24j) is Beer's law. 

To understand what determines the validity of Eq. ( l2T|) . and thus the applicability of 
Beer's law, we calculate, using Eqs. ( ITSll . (fTTl) . and ( fT8l) . the Fourier transform of the 
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Auit 



Ep-Ei-Lo 



+ KEUv(a:, uj + Ei-Ep). (25) 



EUV-induced polarization: 

/oo 
dtP^^\x,t + x/c)e 
-oo 

I J' F 

^ E*p - Er + u ^ 

Therefore, the ratio P^^\x,ijj) / 8-£,\jy{x,ijj) is independent of the EUV electric field — i.e., a 
well-defined linear susceptibility is obtained — only if 

^Euv(a:, UJ + Ei- Ep) 

= = const. (26) 

£^-yjy{x,u:) 

This condition is always satisfied if only terms with I' = I contribute to the polarization 
response in Eq. (125|) . Generally, this is not the case. Still, Eq. ( 126|) serves as a good 
approximation if the EUV pulse is much shorter than 27t /{Ej — Ej/) for /' ^ I. Furthermore, 
for Beer's law to be rigorously applicable, this short pulse must not undergo any significant 
distortion as it propagates through the dense gas, so that the EUV pulse remains short in 
comparison to the dynamical time scales characterizing the target medium prepared by the 
NIR pump pulse. If Eq. (126|) is not satisfied, the propagation equation (!20|) must be solved 
numerically. In the following sections, we discuss both regimes in more detail. 



C. The configuration expansion 

So far, our treatment of the electronic-structure problem has been general. In order to 
describe resonant EUV absorption by laser-generated ions, we now adopt the mean-field 



model discussed in Ref. 



36| and write the NIR-only state vector after the NIR pulse in 



terms of Slater determinants: 

|*NiR,t) = J]a,e-^(^^-^«)*|/) (27) 

+ Z1Z1"*"^~'^'"~''^*I*^)• 
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In this approximation, the set {|/)} contains the ground-state determinant |$o) and all 
particle- hole configurations |$°) obtained from |$o) by exciting or ionizing an electron from 
an occupied spin orbital i (a "hole" orbital) to an unoccupied spin orbital a (a "particle" 



orbital) 37H39|. A particle- hole configuration corresponds to an excited or ionized electron 
plus an ion core with a hole in some shell that is fully occupied in the ground state of the 
neutral atom. The orbital energies of the hole and the particle are denoted by Si and 
respectively. The numerical calculation of the coefficients o^t in Eq. fl7r|) is described 



in Ref. |36|. Since Z is a one-body operator, the sum over the states \F) in Eqs. (IT^ 
and (125|) extends over |$o); the particle- hole configurations, and the two-particle-two- hole 
configurations. The various terms that arise in this way describe processes such as EUV 
absorption by the neutral ground-state atoms and EUV absorption by the ion core. 
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FIG. 1: Atomic levels of Kr'^ associated with resonant absorption at photon energies near 80 eV. 
The notation nlj^ indicates that, relative to the closed-shell ground state of the neutral atom, an 
electron is missing (a hole is present) in the nlj subshell. 



In the following, we focus on the experiment of Ref. js] (see Sec. H]), which studied the 
structures in the EUV absorption spectrum of strong-field-generated Kr"*" ions associated 
with exciting an inner-shell 3d electron in Kr"*" into the outer-valence 4p vacancy created 
by the NIR pulse. Resonant 3d — 4p photoabsorption by the Kr"*" ion core was used to 



10 



study the properties of the hole generated via strong-field ionization. Since such transitions 
involve only hole orbitals, we refer to them as hole-hole transitions. The relevant transitions 
are indicated in Fig. [TJ The transition energies shown in Fig. [T] were calculated using the 

n 

multiconfiguration Dirac-Fock program package GRASP [40']. The configurations that were 
included in the calculation are 4:p~J^^, 4:p~^^, 4:S^^^, 3(igy25 ^.nd Sdg^g- 

We neglect that the EUV field might induce particle-particle transitions. This is an 
excellent approximation, because the NIR-generated photoelectron interacts only weakly 
with the EUV field. Since the electron excited or ionized by the NIR pulse is a spectator in 
the hole-hole transitions of the ion core, one can describe the ions with a reduced density 
matrix 36|: 

pSr^^) = e"^'--'^'^' «' = ^^'-''^'f^f^ (28) 

a 

where the summation is performed over all unoccupied orbitals. Since we consider the 
regime where the NIR and EUV pulses do not overlap, the auxiliary matrix p-)°"^ is time- 
independent. 

The polarization response of the ions to the field of the EUV pulse can be expressed via 
the density matrix by inserting the ansatz (1271) for | \E'nir, t) into Eq. ( lT2l) and following the 
same steps in the derivation that led to Eq. ( l25l) . This yields 

p(^) (X, c) = J] pSr^ J] .,,.4 

>£-E\jv{x,u + Si- Si'). (29) 



Ei - + w 



Here, the absorption of an EUV photon fills a hole in orbital i with an electron from orbital 
z, which has the orbital energy e^. This process is described by the dipole matrix element z-;^. 
The second dipole matrix element, describes the process upon which the hole created 
in orbital % is filled by an electron from orbital i! . 

In the following, we express the hole orbital energies in terms of ionization potentials, 
li = —Ei being defined as the minimum energy required to create an ion with a hole in orbital 
i from a neutral atom in its ground state. Furthermore, since 

p{i) 

[X, —u ) = (x, 00)) , it 

is sufficient to calculate the polarization response only for positive frequencies. Dropping the 
counter-rotating term in Eq. ( l29l) and using ionization potentials instead of orbital energies. 
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we obtain 

P(i)(x, a; > 0) = nAT V pjf^ V ^ 4uv(x, uj-k + I,), (30) 

where T-- is the decay width of the one-hole channel i. We emphasize that Eq. ( l30l) is valid 
only if the EUV pulse comes after the NIR pulse. Recall also that the NIR pulse is assumed 
to remain unmodified as it propagates through the gas. If we dropped this assumption, then 
the ion density matrix would depend explicitly on the atomic position x along the NIR pulse 
propagation axis. The computationally expensive electronic wave-packet problem (I3l) would 
then have to be solved not only once, but for every grid point used in the discretization of 

X. 



III. SHORT-PULSE APPROXIMATION 



We now analyze the EUV-induced polarization assuming that the probe field may be 
approximated by a delta function centered at t = tEuv, i-e., 

SEUv{x,t + x/c) (X S{t -tEUy). (31) 

Within this approximation, the Fourier transform of the EUV electric field is 

SEUY{x,oo)<xe'^'^^^. (32) 
Therefore, we obtain from Eqs. fl2T]) . fl28|) . and fl30|) the linear susceptibility as 



nAT > Pw [tEVY) > 7r. r 



UJ 



2 

This equation shows that the linear susceptibility has poles that are simply related to tran- 
sition energies of the ion core. The existence of an x-independent linear susceptibility allows 
one to use Beer's law [Eq. (Elj)] to calculate the spectrum of the EUV radiation transmitted 
through the target medium. From a practical perspective, this means the following: Under 
the conditions assumed, one obtains, by taking the logarithm of |£^euv(-^) i^)P/|^^euv(0, i^)P, 
a quantity that is proportional to the EUV one-photon cross section 



a«(c.) = 4vr^^^^^!^^ 



(34) 

uat 
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Now we are ready to give an explicit expression for the EUV one-photon cross section of 
strong-field-generated Kr"*" ions. Let j be the total angular-momentum quantum number of 
an orbital hole created by the NIR pulse, and let m be the corresponding projection quantum 



number. We exploit the fact that the reduced ion density matrix is diagonal in m |36| and 
denote the ion density matrix elements by Pj"v • Furthermore, it can be shown that terms 
containing p^- give the same contributions to the polarization response as those containing 
p^j^} ■ We refer to the off-diagonal element P3/2^]/2(^euv) of the density matrix between the 



3/2' 



m 



-1/2 and the ApA 



/2' 



m 



+ 1/2 ionization channels as the coherence. This 



element is equal to p^y2^{/2e~'''''''/' [cf- Eq. ([2SD], i.e., |p^72'v2(W) 
const. The complex constant generally differs from zero, unless the statistical mixture 

described by the ion density matrix is completely incoherent. For a perfectly coherent 



,(1/2) 



|.(l/2) I 
1^3/2,1/21 



hole wave packet, \pf/2i/2\ 



Ref. 



would equal 

n 



^3/2^3/2^1/2^1/2- demonstrated theoretically in 



36| and experimentally in Ref. [8|, strong-field ionization does not in general produce 
•fectly coherent hole wave packets. In other words, generally \p^^f2i/2\ \J ^3/2^3/2/^ 1/2^1/2- 
Using Eq. (IMl) and the notation just introduced, the EUV one-photon cross section 



associated with the hole-hole transitions reads, in the case of krypton. 



47r— Im< 

c 



1(4^3-/2111^113^-/2)1' 



T , A^Zd 



'4P3/2 

(3/2) _2_ 
^3/2,3/2 



1(4^3-/21115113^3-/2)12 



(1/2) 



1 



3/2/ 



'3(i. 



'3/2 



— UJ 



+ 



2 -'4p3/2 

|(4p-/2||D||3(i^ 



(3/2) _3_ 
P3/2,3/2-^Q 



(1/2) 
^3/2,3/2 



30 



3/2)1' 



'3d 



■3/2 



_ „-r3d _ r 

' 2 -'4Pi/2 



.(V2) . 
Pl/2,1/2. 



3\/10 



(4^)372 1 1 D 1 1 3^37) (4^-/2 1 1 D 1 1 3^3-/2) 



''3/2 



3/2/ 



P3/2,l/2l*EUVj 



'3d. 



■3/2 



■4pi/2 



UJ 



+ 



(4/2^1/2 (^EUV))' 



'3d 



'3/2 



'4p3/2 



UJ 



(35) 



Here, (4pJ^| |D| |3(i7) is a reduced dipole matrix element 4l|-|43| (not to be confused with 
a matrix element of the reduced density matrix of the ion). Equation ( 135|) consists of four 
distinct terms. The first three terms are independent of the time delay (because the hole 
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populations pff after the NIR pulse are constant) and describe Lorentzian line shapes 
associated with the three resonant transitions indicated in Fig. [1] The fourth term is a sum 
of absorptive and dispersive line shapes and depends on the coherence P3^^/2(^euv)) which 



is a periodic function of the pump-probe time delay. The period (6 fs 35|) is defined by 
the energy difference between the 4p~^2 ^^^d 4p~^2 channels. Note that for the hole dynamics 
to be observable in the spectrum of the transmitted radiation, it is not necessary for the 
relevant resonance lines (4^3/2 ~^ 3(i~^2 ^^'^ ^Pi}2 ~^ 2*^3/2) spectrally overlap. Only 
the coherent excitation of the two resonances is required, which we achieved by using an 
attosecond probe pulse, i.e., a pulse with sufficient coherent bandwidth. 

Once the EUV one-photon cross section is measured as a function of the photon energy 
u and the time delay tEuv? one can use Eq. ( l35l) to extract all nontrivial entries of the ion 
density matrix [^3/2^3/2; P3/2^3/2p-/'i/2^i/2' ^^"^ P3/2^i/2(^euv)]- This is the basic idea underlying 
the analysis presented in Ref. |8|. Such an analysis allows one to characterize (a) the degree 
of alignment of the j = 3/2 level by comparing pff2l/2 ^^'^ /'3/2^3/2 [^^^ system is fully aligned 
if P3/2^3/2 = ^iid P3/2^3/2 7^ 0) is Unaligned if pff2l/2 ~ ^3/2^3/2]' 0^) population of 

(1/2") (3/2) (l/l) 

the j = 1/2 level relative to the population of the j = 3/2 level [Pi/2 1/2/(^3/2 3/2 ~'~ '^3/2 3/2)]' 

(1/2) |, /(1/2) a/2) 



and (c) the degree of coherence by calculating \p]^l2 1/2 1/ y P3/2 3/2P1/2 1/2- 

In order to be able to determine the ion density matrix elements using Eq. fl35|) . the 
reduced dipole matrix elements must be known. To this end, we proceeded as follows. 



401], we calculated the oscillator strengths for transitions from the Sd, ^ 



Using GRASP „ „ „ 

levels to the 4p^^ levels. In this way, we obtained | (4^3^21 1 1 1 3*^5/2 )P — 0.119 a.u., 
I (4^3/!^ 1 1 1 1 3^3/^2) P = 0.0126 a.u., and | (4p-^2l 1^1 |3f^3/2) P = 0-0695 a.u. The relative ra- 
tios are very close to what would be obtained within the LS coupling scheme. This allowed 
us to employ standard angular-momentum algebra within the LS coupling scheme 



43| to determine the relative sign between (4^3^21 1-^1 1 3c?3/2) (4Pi/2l l-^l |3'^3/2)- Thus, 
(4^-^211^113^^3/2)^^211^11^^3/2) = -0.0296. The natural lifet ime broadening from the 
Auger decay of the 3(iJ^ levels is T^d = 0.00323 a.u. 4^. We mention that, when an- 
alyzing experimental transient absorption data, the spectrometer resolution must also be 
taken into consideration. 

Based on the approach just described, a fit was performed in Ref. js] to determine, directly 
from the experimental attosecond transient absorption data, the reduced density matrix 
elements of strong-field-generated krypton ions. The results are collected in Table HI Also 
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TABLE I: Normalized density matrix elements of strong-field-generated Kr+ ions (a) extracted 



from experiment ^] using Eq. (j35p and (b) calculated using the theory from Ref. 
NIR pulse parameters specified in Ref. [8]. 



36[ assuming the 





(a) experiment 


(b) theory 


2^(3/2) 
^^3/2,3/2 

^^3/2,3/2 

^/^l/2,l/2 


0.23 
0.42 
0.35 


0.05 
0.69 
0.26 


^3/2,1/2 


0.12 


0.13 



shown in Table [T] are the ion density matrix elements calculated using the time-dependent 
multichannel theory developed in Ref. 36|. The NIR pulse parameters assumed in the 
calculation were taken from Ref. jsj. Note that the ion density matrix elements have been 
normalized such that the trace of the ion density matrix equals unity. It is evident from 
Table [T] that experiment and theory give a similar degree of coherence. Also, experiment and 
theory give a similar value for the population of the j = 1/2 level relative to the population 
of the j = 3/2 level. However, the experimental data suggest a much smaller degree of 
alignment of the j = 3/2 level than predicted by our theory. 

The attosecond transient absorption cross sections [Eq. (l35l) ] for the two sets of density 
matrix elements are plotted in Fig. [2j Again, there is no overlap between the NIR and 
EUV pulses, that is, the NIR pulse is centered at a large negative value of tEuv- The 
strongest of the three absorption lines does not depend on the time delay. This absorption 
line corresponds to a transition to the 'id'^^^ level. Since this level can be reached only 
from the 4p~^2 level, but not from the 4p~^2 level (see Fig. [1]), the strongest absorption line 
is insensitive to the coherence between 4:p~^^ and 4p~^2- The other two transitions, which 
both involve the 3d~J^^ level (see Fig. [1]), display a conspicuous dynamical behavior as a 
function of the time delay. The modulation of the weakest line, which corresponds to the 
transition 'ip'^J^^ — 3^3^25 is most pronounced due to the coupling to the relatively strong 
4p~y2 3(^3^2 transition. Moreover, for certain delays the absorption cross section for that 
line even becomes negative. This happens when the coherent population transfer from 4p~^2 
to 4p~y2 '^^3/^2 dominates over the absorption from the 4p~^2 state. 

Apart from the oscillation of the respective line strengths with the spin-orbit period of 6 
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FIG. 2: (Color online) Attosecond transient absorption cross section, in Mb, of strong-field- 
generated Kr"*", plotted as a function of the photon energy and the time delay. The cross section 
was calculated using Eq. (j35p . Panels (a) and (b) show the attosecond transient absorption cross 
section for the theoretical ion density matrix elements in Table HI In panels (c) and (d), the attosec- 
ond transient absorption cross section is shown for the experimental ion density matrix elements 
in Table HI 

fs, the energetic positions of the resonance hnes oscillate as well. This is particularly easy 
to see in panels (b) and (d) in Fig. [2j These energy oscillations are a consequence of the 
interplay between the absorptive and dispersive terms mentioned earlier in connection with 
Eq. (ESD. 
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IV. THE ACCURACY OF BEER'S LAW 



In the previous section, our discussion was based on the approximation that differ- 
ent frequency components of the EUV pulse propagate independently. The assumption 

p(i) 

{x,u}) oc SEvy{x,uj) allowed us to describe the polarization response with a linear sus- 
ceptibility x^^H^) and to integrate Eq. (!20!) analytically. In this section, we investigate the 
accuracy of this approximation. 

As follows from Eq. (!30|) . the polarization response P^^\x,u) is not proportional to 
S-evy{x,u}) if the density matrix contains non-zero off-diagonal elements that correspond to 
states coupled by dipole transitions through an intermediate excited state. Physically, this 
means that an ion that absorbs a photon with an energy Ui can coherently emit a photon 
with a different energy U2, provided that the initial ionic state is a coherent superposition 
of two or more states. In this typical A-scheme, the polarization response at the photon 
energy U2 obviously depends not only on £^EUv(a;, 1^2), but also on £-evy{x,ui). In this case, 
Eq. (!2T]) is an approximation. 

To go beyond this approximation, we numerically solve the first-order propagation equa- 
tion ( 12U]) using Eq. ( 15U]) to evaluate the polarization response of the medium at each prop- 
agation step. Even though P^^\x,uj) is no longer proportional to SEvvix.co), we compare 
the results of numerical propagation with those obtained in the previous section in terms of 
the apparent one-photon cross section: 

2 

^. (36) 



uatL 



In Fig. 131 we show a false-color representation of (Tapp(i^) evaluated for different delays 
between the NIR pump and EUV probe pulses. For this simulation, we used a bandwidth- 
limited Gaussian EUV pulse with a central photon energy of 80.8 eV and a full width at 
half maximum of intensity equal to 150 attoseconds. The elements of the density matrix 
were taken from Table [H The propagation in a gas of strong-field-generated Kr^ ions with 
an atomic number density of tt-at = 2.2 x 10^^ cm^^ was terminated after L = 1 mm. The 
spectrum of the EUV pulse before and after propagation, for the experimental ion density 
matrix elements in Table [H is shown in Fig. |H 

A careful inspection of Figs. O El and H] reveals that the approximation underlying Beer's 
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FIG. 3: (Color online) Apparent attosecond transient absorption cross section, in Mb, of strong- 
field-generated Kr+, plotted as a function of the photon energy and the time delay. The cross 
section was calculated by applying Beer's law to the numerically propagated EUV field. Panels 
(a) and (b) show the apparent attosecond transient absorption cross section for the theoretical 
ion density matrix elements in Table HI In panels (c) and (d), the apparent attosecond transient 
absorption cross section is shown for the experimental ion density matrix elements in Table HI 

law notably affects the line at 81.1 eV, which corresponds to the transition 4p~^2 ~^ "^'^3/^2- 
This is not surprising: for the same reasons why this absorption line exhibits strong quantum 
beats, it is also sensitive to other effects related to off-diagonal elements of the density 
matrix. Still, in spite of the strong absorption, the discrepancy between the model assuming 
the validity of Beer's law and the results obtained by numerically solving the propagation 
equation is rather small. 

It is instructive to repeat the analysis based on Eq. (!35|) and retrieve the density matrix 
from the apparent absorption cross section as if Beer's law were rigorously valid. Table Hi] 
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FIG. 4: (Color online) The spectrum of the EUV pulse before (dash-dotted line) and after prop- 
agation for tEUV = 0. The dashed line, representing a simulation assuming the validity of Beer's 
law, corresponds to panels (c) and (d) in Fig. [2l The solid line shows the result of numerical 
propagation, corresponding to panels (c) and (d) in Fig. [3l 



TABLE II: Apparent density matrix elements of strong-field-generated Kr"*" ions extracted from 
the apparent attosecond transient absorption cross section for (a) the experimental ion density 
matrix elements in Table [I] and (b) the theoretical ion density matrix elements in Table HI 





(a) 


(b) 


2„(3/2) 

^^3/2,3/2 
^^^1/2, 1/2 


0.21 
0.44 
0.35 


0.03 
0.71 
0.26 


1^^3/2,1/21 


0.14 


0.14 



gives these apparent density matrix elements. The good agreement between Tables [T] and [TTl 
indicates that Beer's law is indeed a good approximation for extracting electronic-structure 
information from transient absorption data, even though it should be used with care. We 
may conclude, in particular, that the discrepancy between the experimental and theoretical 
degrees of ion alignment (j = 3/2 level) cannot be explained by a failure of Beer's law. 
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V. CONCLUSIONS 



In this paper, we discussed the theory underlying attosecond transient absorption spec- 
troscopy of strong-field-generated ions. This theory was employed in Ref. jsl to experimen- 
tally determine the reduced density matrix of Kr+ ions produced by an intense NIR pulse. 
Good agreement between experiment and theory was found for the degree of coherence and 
for the population of the j = 1/2 level relative to the population of the j = 3/2 level. 
However, experiment suggests strongly suppressed alignment of the j = 3/2 level, which is 



36|. The origin of 



not consistent with calculations based on the theory described in Ref. 
this discrepancy is currently unknown. 

As demonstrated in this paper, EUV propagation effects beyond Beer's law do not explain 
the discrepancy between experiment and theory found for few-cycle NIR pulses jsj. Earlier 
measurements on Kr+ ions generated using 50 fs NIR pulses gave a degree of alignment in 



rather good agreement with an adiabatic strong-field-ionization theory 



16| . A noticeable 



reduction of ion alignment, in comparison to the adiabatic strong-field-ionization theory. 



was observed in Xe"*" ions generated using 45 fs NIR pulses [29|]. But the effect was not as 
pronounced as it is in the current case, and it was surmised at the time that the discrepancy is 
a consequence of nonadiabatic effects 29|]. However, nonadiabatic effects cannot explain the 
disagreement between the experimental and theoretical ion density matrix elements shown 
in Table [U for the strong-field-ionization theory employed 



is based on numerical wave- 



packet propagation and does not suffer from the limitations of the adiabatic approximation. 

It seems likely that the origin of the discrepancy must be sought in either of the following 
two possibilities. The first possibility is that in the experiment of Ref. js], multielectron 
effects beyond the multichannel theory of Ref. 
of Ref. 



played an important role. Since, in view 



16j, these multielectron effects would appear to have a smaller impact when using 



ectron 
8| was 



longer NIR pulses, the observed discrepancy might suggest an enhancement of multie 
effects by few-cycle pulses. The second possibility is that the experiment of Ref. 
affected by substantial NIR propagation effects. This could be clarified by repeating the 
experiment at a lower target density. 
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Equation (21) in Ref. 



cr{u,to) = ^— Re 
he 



22( 1 should read 

dt e-^2itgi^t X (i(to)|e^^i*/V2ie"^'^'*/V2iK(to)> 
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